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Abstract 

Feynman's path integral is herein generalized to the nonextensive canonical 
density matrix based on Tsallis entropy. This generalization is done in two 
ways by using unnormalized and normalized constraints. Firstly, we consider 
the path integral formulation with unnormalized constraints, and this gen- 
eralization is worked out through two different ways, which are shown to be 
equivalent. These formulations with unnormalized constraints are solutions 
to two generalized Bloch equations proposed in this work. The first form of 
the generalized Bloch equation is linear, but with a temperature-dependent 
effective Hamiltonian; the second form is nonlinear and resembles the anoma- 
lous correlated diffusion equation (porous medium equation). Furthermore, 
we can extend these results to the prescription of field theory using integral 
representations. The second development is dedicated to analyzing the path 
integral formulation with normalized constraints. To illustrate the methods 
introduced here, we analyze the free particle case and a non- interacting scalar 
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field. The results herein obtained are expected to be useful in the discussion 

of generic nonextensive contexts. 
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I. INTRODUCTION 



Path integral techniques constitute a powerful tool, employed to analyze countless phys- 
ical situations, for instance, in standard statistical mechanics (Boltzmann-Gibbs). On the 
other hand, there is nonextensive behavior common in many branches of physics, for exam- 
ple, systems with long-range (gravitational) interactions, long-time memory, and fractally 
structured space-time. This nonextensive behavior indicates that the standard statistical 
mechanics and thermodynamics need some enlargement. Consequently, the path integral 
employed in standard statistical mechanics also needs some generalizations in order to ac- 
commodate the nonextensive effects. 

Recently, a theoretical tool based on nonextensive entropy (Tsallis entropy) Q has suc- 
cessfully been applied in many situations, for example, Levy-type anomalous superdiffusion 
@, Euler turbulence ||, self-gravitating systems @], cosmic background radiation [||, pecu- 
liar velocities of galaxies 0, anomalous relaxation through electron-phonon interaction 
ferroffuid-like systems ||, and nonlinear dissipative dynamical systems ||. In this context, it 
is very important to understand the main properties of the generalized statistical mechanics 
based on Tsallis entropy more deeply. In particular, it is interesting to develop basic tools 
for calculations such as the semi-classical approximation [|I(J , perturbation |jTiyi2| and vari- 
ational methods [|TT|-|TE[, linear response theory Green functions |T3|,|TB[, non-Gaussian 
integrals ||17|| , and path integral, which is the subject of this work. 

We develop two path integral formalisms and Bloch equations based on Tsallis statistics 
with unnormalized constraints in Section II. The first form of the generalized Bloch equation 
is linear, but with an effective Hamiltonian, which has a temperature dependence, and 
the other equation resembles the anomalous diffusion equation (porous medium equation). 
The path integral solutions to these two different but equivalent generalized Bloch equations 
are also obtained in Sections II as well as the extension to the field theory. The path 
integral is first formulated for effective Hamiltonian case. In the other case, we use integral 
representations to formulate the path integral. To illustrate theses methods, the free particle 
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case for q > 1 and q < 1 and the free scalar field bosonic are studied. Section III is 
dedicated to analyzing the path integral formulation obtained with normalized constraints 
and its extension to quantum field theory, where an example is exhibited. Conclusions are 
presented in Section IV. 



II. PATH INTEGRAL BASED IN UNNORMALIZED CONSTRAINTS 

Let us start this section with some remarks about Bloch equation and nonextensive Tsal- 
lis statistics in order to introduce the generalized Bloch equations. After this presentation, 
we obtain the path integral formulation for the solutions of these Bloch equations, giving 
as example the free particle case. Furthermore, we apply these developments for a brief 
discussion about quantum scalar fields. 



A. BLOCH EQUATION 

The usual density matrix 

frGS) = e-^/Zi (Trp!(/?) = 1) (1) 

{Z\ = Tr e~^ H is the usual partition function) can be obtained directly from a maximum 
entropy principle |18 , and the unnormalized matrix 



fk(P) = exp(-/3H) (2) 



obeys the Bloch equation |B| (see also Ref. |2(J p. 48) 

(3) 

with the initial condition pi{(3 — 0) = 1. When is calculated and subsequently Z\, we 

can obtain the usual free energy F\ = — (1//3) InZi. 

Recently, the density matrix ([I]) was generalized in the nonextensive context []I|,|21,2"2"|, 



p q ((3) = [l-(l-q)(3H] 1 ^/Z q , (4) 

where Z q = Tr[l — (1 — q)(3H} l ^ 1 ~ q " > is the generalized partition function. In Eq.([|) we are 
supposing that 1 — {l — q)f3E n > 0, where E n are the eingenvalues of H. When this condition 
is not satisfied, we employ p(E n ) = 0, i.e., a cut-off is used so that p q (E n ) remains positive. 
By using Z q , we can obtain the generalized free energy F q = — {l/(3)(Z q -~ q — 1)/(1 — q). 
These generalizations were obtained by employing the maximum entropy principle with the 
Tsallis entropy, 

S q = -k Tr tl£. , (5) 
1 — q 

and the constraints (H) q = Tr p q q H and Tr p q = 1 (k is a positive constant). The parameter 
q G R gives the degree of nonextensivity, and in the limit q — > 1 the usual statistical 
mechanics is recovered. Note that the constraint of the energy is unnormalized, since (l) q = 
Tip q 7^ 1. We return to this point at Section IV (see Ref. [^J for a complete discussion 
about the choice of the constraints). 

Following the previous discussion about the Block equation, we employ the unnormalized 
matrix 

Pq ((3) = [l-(l-q)PH] 1/ ^ ) . (6) 



This matrix satisfies alternatively the equations 

dp q H 



p q , (7) 



OP 1 - (1 - q)pH 

- 8 -k = Hph (8) 



with the initial condition p q (j3 = 0) = 1. These equations generalize the Bloch equation 
in the nonextensive Tsallis context. In fact, the generalized Bloch equation (|7|) is linear in 
p q and has the same form of the Eq. (|3|), if we employ the effective Hamiltonian H e ff = 
H[l — (1 — g)/3i/] _1 . On the other hand, the generalized Bloch equation (^) has some formal 
resemblace to the anomalous correlated diffusion equation EM. In the following, we develop 
the path integral solution for both generalized Bloch equations. 



B. PATH INTEGRAL FORMULATION WITH H eff 



Since Eqs. (|3]) and (|7|) have the same form, the path integral representation for p q has the 
same structure of the usual path integral. Thus, the path integral for p q can be written as 



p q (x,x';(3) = (x\p q ((3)\x') 



lim 

N^oo 



T>xT>p exp 



7V-1 

n dx 

n=l 

1 rP 



i 2nh 

±J P (ipx-H e „)dr 



(1 N 

ex P |^ Pn{ x n ~ Zfi-i) - eH eff (p n ,x n ,ne)} 



(9) 



where x n = (x n+ i— x n )/2, x = xo, x' = xn and e = /3/iV. In this path integral representation, 
we are considering that the Hamiltonian H as well as H e f / are Weyl ordered (if this is not the 
case, the discussion has to be redone along the present lines), and consequently the midpoint 
prescription was employed |25 |. In general, this path integral is not easy to calculate because 
the Hamiltonian H e ff(— H[l — (1 — q)^!!]^ 1 ) has an unusual nonlinear dependence on the 
moments and coordinates. However, as we shall see, it is possible to calculate p q for the 
free particle case. Simple textbook examples, such as free particle, the harmonic oscillator, 
the non- interacting scalar field, etc., do not exibit any kind of nonextensive behavior. It is 
not necessary to introduce a generalized statistical mechanics to deal with these systems. 
However, these simple examples are very useful in order to illustrate how the formalism 
works. And they are even more instructive when they provide exactly solvable cases. In this 
case, considering the free particle case, we first perform the integration in the coordinates, 



p q (x,x';(3) 
x exp <^ -e 



dp 
2nh 



exp 



(x — x 



lim 

N^oo 



N-l 



Jl dp n 5(p n -p n . 



n=l 



N 



n=l 



p 2 n /2m 



1 — en(l — q)Pn/2rn 



(10) 



By using the identity (see appendix A) 

ln(l + x)= lim V ( ^ ) —^—r 

V ; 1 \NJ 1 + f 



N 



(11) 



and considering the case q > 1 with (3 > 0, we verified that Eq. (|T0|) reduces to 
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p q (x,x';(3) 



roc dp 


'ip f a" 




/ o % ex P 


—{x - X ) 





1 - (1 - q)(3 



p 

2m 



1/(1-3) 



2m 



ir(q-l)phr 



1/2 



m(x — a;') 2 \ ^ 2 
2{q-l)ph 2 ) 



i/(g-i)-l/2 



x 



'2m{x — x') 2 \ 1 ^ 2 



(12) 



( g -i)^ 2 ; 

In the last step, we used a known integral (for instance, see Ref. [2~B| p. 321) and relation 
between the Wittaker function Wq^ v and the modified Bessel function of second kind K v . 

In the case q < 1 with (3 > 0, that has a cut-off, the integration limits of the moments 
are limited and the path integral (|10D reduces to the first part of Eq. (|l^) with p satisfying 
the condition 1 — (1 — q)(3p 2 /2m > 0. Thus, when the remaining integration is performed 
(see Ref. |26|] p. 321), we obtain 



p g (x,x';(3) 



■ 2m I 1 / 2 j 

(T^pJ dp 



i 1/2 2%h 
+ 1 



exp 



. (.1-1)0 
' 1 



— [x — X , 
n 



1 



1/(1-9) 



\l-q 

X Jl/(l-q) + l/2 



m 



nl/2 



2tt(1 -q)(3h 2 _ 
2m(x — x') 2 \ 



' 2(l-q)Ph 
^ m(x — x') 2 



2\ 1/2' 



l/(l-g)+l/2 



X 



(13) 



1 - g)/^ 2 

The partition function can be obtained for q > 1 and q < 1 by taking the trace of 
p g (x,x'; (3). In fact, for g > 1 we obtain 

rL/2 



L/2 



dx p q (x, x; (3) 



m 



lV2 tv_J I' 

Vo-l 2- 



27r(g-l)^ 2 J r(-ij) 



(14) 



and for g < 1 we have 



Z 9 = L 



??? 



i 1/2 ni+i; 



(15) 



.27r(l-g)^ 2 J T^ + f) ' 

where L is the length of the integration region (box size). Because L is large the expressions 
( PI) and ( |T5D reduce to the generalized classical partition function. 
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C. PATH INTEGRAL FORMULATION WITH INTEGRAL REPRESENTATIONS 



Let us consider the Hilhorst's identity (private communication to Tsallis 

[1 - (1 - q)pH} 1,{l ~ q) = —, v / dv v—^e-v exp(-v(q - 1)/3H) . (16) 

r (ft) Jo 

This operator identity is essentially a direct application of the usual integral representation 
of the gamma function. Indeed, if both sides of the Eq. ([16]) are applied on an eigenstate 
of H with q > 1, (3 > 0, and E n > (E n is an arbitrary eigenvalue of H) the identity ( |T6"D 
reduces to the usual integral representation of the gamma function. From Eq. ([IB] ) the path 
integral representation for p q (f3) can be obtained. In fact, it is necessary to employ only the 
path integral representation for f>\{v(q — l)/3) in the last part of Eq. ([16|). Thus, 

p q (x,x;(3) = , - r / dvvi- 1 e~ v pi(x,x ;v (q - 1)(3) , (17) 

where pi(x,x';v(q — l)f3) must be replaced by its usual (or holomorphic) path integral 
representation. Furthermore, if we take the trace of Eq. (0), we obtain the connection 
between the generalized and usual partition functions P7| . 

As illustration , let us apply Eq.(|T7D to obtain the generalized density matrix for a free 
particle in the one-dimensional case. For this example, the usual density matrix in the 
coordinate representation is (for instance, see Ref. [^0| p. 49) 



M^/W^exp 



m 



i\2 



(x — x') 



(18) 



When this expression is replaced in Eq. ( |T7D and the integral is calculated (see Ref. [EH] p. 



935), we obtain, as expected, Eq. ( |12D 



The integral transformation ( |16|) can not be applied for q < 1, because the integral is 
not convergent. However, for q < 1 the Prato's identity 

[1 - (1 - q)f3H} l/{1 - q) = — r [ ' 2 —^\ [ dv (-vy^e- v exp(v(l - q)(3H) (19) 

27r \1 — q J Jc 

can be employed []2S|] , except for some particular values of q, (2 — q)/(l — q) = integer. As 
in the case q > 1, this integral representation is a direct consequence of another gamma 
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function representation (for instance, see Ref. p6[ p. 935). The contour integration can be 
chosen in order to eliminate the condition (2 — q)/(l — q) 7^ integer, see for instance, Refs. 
fT2|| and | 2"9]| . Thus, we can formulate a path integral for q < 1 using Eq.(|T9"D. In fact, it 



is sufficient to replace exp(t> (1 — q)(3H) by its path integral representation in Eq. fllTf ) to 
obtain the desirable path integral representation for p g (x,x'; (3) with q < 1. 

D. EXTENSION TO THE FIELD THEORY FORMALISM 

We can use the previous path integral prescription to study some aspects of the Tsallis' 
nonextensive thermostatistics as applied to field theory. Let us start our discussion about 
nonextesive statistical field theory by considering the scalar field. In this example, we obtain 
the generalized partition function and the pressure attributed to the ground state of a free 
scalar field. The Lagrangian density is 

C = \d^d^ ~ \m 2 <P , (20) 

with c = 1, h = 1, and consequently the Hamiltonian density is 

H = + \ (V0) 2 + . (21) 

By taking the trace of the Eq.(|H5|), the generalized partition function for the free field, we 
get 

Z q = j c dv Kf\v) J---J VcffDix ef ^Jd 3 x (^- ; (22) 

with (3* = (1 - q)(-v)P and = (f=j) (-vy(^e~ v . This partition function, after 
the momentum integration reduces to 

Z q = J dv Kf\v)f---J V<\> efo* d - / d 3 ^^*) . (23) 



By substituting the Lagrangian density and following the usual calculation [j3*0[| , we obtain 
for the generalized partition function the following expression 
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(24) 



where u = (p 2 + m 2 ) 1 ^ 2 and uj n = 27m/ P*. Performing the sum over n, the above expression 
can be rewritten as 



Z q = J o dvKW{v)e*pl-vJ ^[^{l-q){-v)u + \n(l-e-^^) 



(27T) 



(25) 



By taking m — > this expression recovers essentially the generalized partition function for 
the blackbody radiation [ 29] . In the following, as illustration, we calculate the generalized 
pressure attributed to the ground state. In this case, we use the definition of the generalized 



pressure [p9|pT 



Pn 



ld_ ( Z\~*-Y 
(3dV\ q-1 , 



(26) 



In this way, by using Eqs. (p5| ) and ( 2~6|) with T — 1/(3 sufficiently small (a detail discussion 



of this limit is presented in Ref. |12|), we obtain 



X 



X 6XP I ~ V I ^ " q){ ~ v)u ° + ln ( X - e-M-M-***) 



(2tt)' 



d 3 p 



(27) 



where 



Kf\v) 



2ni 



(28) 



At this point, we remark that the result for the generalized pressure of the ground state 
coincides with the same q = 1 case [EtJ. In general, the ground state properties can be 
obtained by using the small temperature limit of the generalized statistics independently of 
the q value. 
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III. PATH INTEGRAL BASED ON NORMALIZED CONSTRAINTS 



As in the above discussion, we can formulate the path integral for nonextensive systems 
using another version of the generalized statistics with different constraints(a detailed dis- 
cussion of the choice of the constraints is presented in Ref. ||23|| ). Considering the case with 



the normalized constraint ((H)) q = Trp H /Trp 9 , instead of (H) q = Trp H, we obtain that 



Pq 



Pa 



(29) 



with p q = (l — (1 — q)ft(H — ((H)) g)^ ^ Q \ Z q = Trp q , and ft = ft/Tip 9 . In the same way 
of the unnormalized constraint case, we can express Eq.(^) through the path integral, 



p q (cf>,<f>',ft) = / duK^i 



c 



T><f)Vir exp 



drd 3 x(<i)-n -H(<f),n)) 



(30) 



where ft = (1 - q)ft(-v) and K^(v) = K^^e^ 1 -^^" . 

To exemplify the normalized formalism, let us consider the free particle case, 

dr[xp — H(x,p)) 



f a dvKV>(y) J ■ 


■ ■ J T>xV>p exp 


[1 









(31) 



Here, the calculation resembles the unnormalized case. Thus, we restrict ourselves to pre- 
senting the main results. For instance, in the q < 1 case, after some calculations, the density 
matrix for one particle is given by 



Pq(x,x',ft) 



mTrp 9 


1/2 " 


_2tt(1 - q)fth 2 _ 





x • I -I Jl/{l-q)+l/2 



1-q 

and, for the q > 1, it reduces to 



mTrp (x — x') 2 v 

2m(x — x') 2 Trp 9 
(1 - q)fth 2 



Tip 9 



—+- 

l-q T 2 



3 

Tip 9 



1/2- 



((H)), 



(32) 



Pq(x,X,ft) 



3-8 
1-9 



ir(q-l)ptf 



X r ^^ l/(g-l)-l/2 



Tip 

2m(x — x') 2 Trp 9 
, (q-lW 



A2m ~9N 1/2- 



/ m(x — x') Trp 

{ 2( q -i)fth 2 

(1 ~q)P. 



_L 1 

q— 1 2 



Tip 9 



■((H)), 



(33) 
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Here, Z q and ((H)) q for q < 1 are 



I 2mw \ 1 ' 2 r(j^) a An(i-,w/; 

£ U -«)/»' J r(f^ + il (1 + (1 ?) 2 ) 



2/(1+9) 



«#»9 = , 

and, for the case q > 1, are 



1 



/ 2m7r \ 1/2 r - |) i 1/(1 _ ff)+1/a 

r(^) 1+1 gj 2 j 



2/(1+9) 



((^))9 = ^Trp g 



In these calculations the general relation Tip 9 = Z\ q was employed. 



(34) 



(35) 



A useful guide to calculate thermodynamics quantities in the normalized formulation is 



based on the relation 23 



if 



Tr~p q q + (l-q)fl™((H)) q ' 
where the f3 unnor aiid /3 nor are Lagrange multipliers associated with unnormalized and nor 
malized constraints, respectively. 



(36) 



IV. CONCLUSIONS 

In this work, the path integral formulation for nonextensive systems based on Tsallis 
entropy was developed, considering both the unnormalized and the normalized constraints. 
For unnormalized formalism, we deduced, in a unified way, two different but equivalent 
generalizations of the path integral as well as their related Bloch equations. In the first case, 
the Bloch equation has the merit of preserving the form when compared with the usual Bloch 
equation, but containing a temperature-dependent effective Hamiltonian. This temperature 
dependence makes the evaluation of the respective path integral more difficult. The other 
generalized Bloch equation has a usual form in the sense that there is no explicit temperature 
dependence; but on the other hand, it is a nonlinear operator equation (it resembles the 
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anomalous correlated diffusion equation). The corresponding path integral (one and the 
same from both generalized Bloch equations) is obtained from integral representations of 
the gamma function. Thus, the usual path integral (q = 1) is directly used to perform its 
generalization (q ^ 1). This fact provides a natural way to generalize the main properties of 
the usual path integral, covering from the usual quantum mechanics to quantum field theory. 
Consequently, the second path integral formulation, based on integral representations, is 
more interesting than the first one. In addition, we analyzed the path integral formulation 
using the normalized constraint in the same way as the unnormalized one, where the free 



particle case was analyzed and our result are in agreement with the classical context p2| for 
one particle in one- dimensional space. Furthermore, the free scalar field was considered. 
In particular, this work can be useful in the future investigations of macroscopic proper- 



ties of one-electron states systems with a multifractal internal structure j33| . This idea must 
be explored in another opportunity. Let us conclude by saying that we expect the above 
generalizations to be useful in the discussion of generic nonextensive systems, hopefully in a 
similar way Feynman's path integral is very useful in many branches of physics. Moreover, 
the present generalization appears as a kind of extension of the scale-invariant diffusion in 
space-time (see figure of [^t| p. 177), in the sense that it could be anomalous (q ^ 1), 
instead of the usual Brownian one (q = 1). 
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APPENDIX A: 

Since we have not found the identity 

ln(l + x)=Jim (*>-!) (Al) 
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in anywhere else, we establish it in the following. In order to retain the convergence in all 
the steps of the subsequent calculation, we employ a new variable y, defined by the relation 
x — m — 1 + y, where m > 1/2, and | y \< m. Thus, the identity (|Al|) can be written as 



N 



ln(m + y) = lim 

N— >oo — ' 



n=l 



m — 1 + y 
N 



1 



n(m— 1+y) 
TV 



(A2) 



After a simple algebra, we obtain 



N 



lim 



n=l 



m — 1 + y 

mN 



kk 



E(-lra 



(A3) 



k=0 



for r.h.s. of 



where a = + ^ n J^- < 1. Furthermore, by using a binomial 



m N m N 

expansion for a, in the above expression we have 



oo k 



N^oo 



fc=0s=0 



m — 1 + y 

mN 



y s (m — 1) 
(mN) k 



N 



Y,n k ~ s {n-N) s . 

. n=l 



(A4) 



To evaluate the sum in n, it is convenient to employ a further binomial expansion, 

, (A 



In 



u=0 



n u (-N) s - u . 



(A5) 



After this, we can sum over n, retaining only the dominating contribution in N, as the other 
contributions vanish in the limit iV — > oo. In this way, r.h.s. of is reduced to 



OO 1 

E 



k=0 



k+1 



m — 1 + y 



m 



k+l 



£(-l) s (m - 1) s y x ' 

s=0 



(A6) 



Performing the sum over s = 0, we get for this expression 



E 

k=0 



■l) k fy\ k+1 ~ 1 /m-l^ fc 

+ E 



fc + 1 Vm 



k=l 



k V m 



(A7) 



The resultant sums are respectively the series representations for In + and ln(m), 
establishing, therefore, the identity (|A2|) . 
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